Adiabatic quantum pump in the presence of external ac voltages 
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We investigate a quantum pump which in addition to its dynamic pump parameters is subject to 
osciUating external potentials applied to the contacts of the sample. Of interest is the rectification 
of the ac currents flowing through the mesoscopic scatterer and their interplay with the quantum 
pump effect. We calculate the adiabatic dc current arising under the simultaneous action of both 
the quantum pump effect and classical rectification. In addition to two known terms we find a third 
novel contribution which arises from the interference of the ac currents generated by the external 
potentials and the ac currents generated by the pump. The interference contribution renormalizes 
both the quantum pump effect and the ac rectification effect. Analysis of this interference effect 
requires a calculation of the Floquet scattering matrix beyond the adiabatic approximation based 
on the frozen scattering matrix alone. The results permit us to find the instantaneous current. In 
addition to the current generated by the oscillating potentials, and the ac current due to the variation 
of the charge of the frozen scatterer, there is a third contribution which represents the ac currents 
generated by an oscillating scatterer. We argue that the resulting pump effect can be viewed as a 
quantum rectification of the instantaneous ac currents generated by the oscillating scatterer. These 
instantaneous currents are an intrinsic property of a nonstationary scattering process. 

PACS numbers: 72.10.-d, 73.23.-b, 73.40.Ei 



I. INTRODUCTION 

Dynamical transport in mesoscopic structures attracts 
presently considerable attention—"—-. In particular the pos- 
sibility to vary several parameters at the same frequency 
but different phases- of a mesoscopic system opens up new 
prospects for the investigation of quantum transport. Apply- 
ing two slowly oscillating potentials at frequency oj with fixed 
phase lag A(p to a mesoscopic conductor connected to reser- 
voirs having equal electrochemical potentials one can generate 
an adiabatic dc current 

/dc ~ <.^sin(A(/3). (1) 

Such a current was measured experimentally-. However 
the precise origin of the measured current is still unclear. 
At least two mechanisms considered in the literature can 
contribute to the experimentally measured current. First, 
there exists a quantum pump ejffeciS*— ~— - which is due to 
quantum-mechanical interference and dynamical breaking of 
time-reversal invariance. Second, there also exists a rectifica- 
tion of ac currents by the oscillating scattere»— if it is 
part of an external circuit with non-zero impedance. Closely 
related to this second effect is a pump in the presence of in- 
elastic scattering: in addition to the externally driven pump 
parameters, inelastic scattering leads to an effective oscillat- 
ing (electro-)chemical potential of the pump which acts like 
an additional pump parametei42i. 

The aim of the present paper is to investigate both these 
mechanisms on the same footing. To this end we consider a 
phase coherent oscillating scatterer coupled to reservoirs with 
oscillating potentials (see Fig. 1). We will show that in gen- 
eral the above mentioned mechanisms do not simply add but 



interfere between themselves. This leads to a renormalization 
of both the quantum pump effect as well as the rectification 
effect in the total dc current. 

Theoretically quantum pumps have been investigated 
mostly under the (implicit) condition that the external circuit 
exhibits zero-impedance. The work of Brouwer'*^, Polianski 




FIG. f : A mesoscopic pump with scattering matrix S(t) os- 
cillating with frequency u is coupled to Nr reservoirs with 
electrochemical potentials ^a{t) oscillating with the same fre- 
quency UJ. A quantum pump effect and a classical rectification 
effect together result in dc currents la flowing through the 
scatterer. The full current is time-dependent and is needed 
to characterize pumps in a non-zero impedance external cir- 
cuit. 
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and Brouweu— , the work on inelastic scattering— mentioned 
already and the recent work of Martinez-Mares, Lewenkopf 
and Mucciolo— represent the few exceptions. In reality the 
zero-impedance condition seems never exactly fulfilled. Cou- 
pling an oscillating gate voltage to a scatterer leads, due to 
the long range nature of the Coulomb interaction, effectively 
to oscillating voltages at all terminalsi?. In addition, in ex- 
periments, the pump is investigated with an impedance in 
series with the oscillating scatterer. Furthermore, a voltage 
probe, to maintain zero current in the presence of pumping, 
in effect generates an oscillating potential which acts back on 
the pump--. Therefore, for theory to make contact with ex- 
periment, it is necessary to consider the effect of oscillating 
voltages at the contacts of the conductor. 

The paper is organized as follows. In Sec^we develop the 
Floquet scattering matrix approach for ac quantum trans- 
port through the nonstationary (oscillating) scatterer in the 
presence of oscillating reservoir potentials. A full theory re- 
quires even to first order in frequency an investigation of non- 
adiabatic corrections to the adiabatic (frozen) scattering ma- 
trix. These corrections are discussed in Sec lIIII The cur- 
rent to linear order in the reservoir potentials is calculated 
in Sec If V I We illustrate the results for a simple one-channel 
scatterer with two contacts. In Sec[3 we present a general 
expression for the current valid for finite potentials. Sec. I VII 
gives the expression for the instantaneous current. 



II. GENERAL APPROACH 

For simplicity we consider here the mesoscopic sample, the 
pump, connected to Nr reservoirs via single channel leads 
FigQ We are interested in the dc and ac currents flowing in 
the system if this system is subject to a cyclic evolution with 
period T. The general situation we want to consider admits 
the scatterer and the reservoir properties to be oscillating with 
frequency uj = 2n /T . 

We use the scattering matrix approach to ac transport in 
phase coherent mesoscopic systems-. According to this ap- 
proach the currents flowing in the system are determined by 
the scattering of electrons coming from the reservoirs by the 
mesoscopic sample—*--. In the present paper we deal with 
non-interacting electrons. A full theory has eventually to 
treat the internal potential of the pump in a self-consistent 
manner. 

The scattering properties of a mesoscopic sample oscillating 
with frequency uj can be described via the Floquet scattering 
matrix S'p (see, e.g., Ref. Isoh . 

The matrix element Sp^^fj{En, E) is a quantum mechanical 
amplitude for an electron with energy E entering the scatterer 
through lead /3 to leave the scatterer through lead a having 
energy En = E + nhui. We use Greek letters a, P to number 
the leads connecting the scatterer to the reservoirs: a,l3 = 

I,..., Nr. 

Denoting by a an annihilation operator for incoming par- 
ticles we can write down the expression for the annihilation 
operators o for outgoing particles— 



(but not necessary stationary) state. Therefore the proper- 
ties of incoming particles are independent of the scatterer and 
are determined by the reservoirs. To be definite we suppose 
that the cyclic evolution of any reservoir a is due to solely an 
oscillating electrochemical potential ^dt)- 



lla{t) = HO, a +eVa{t), 
Va{t) = Va COs{ljjt + <Pa), 
eVa < /iO,Q- 



(3) 



We emphasize that we must keep track of the phase shifts 
ifa since there are a number of oscillating potentials and we 
can not eliminate all the phases ipa simultaneously by merely 
shifting the time origin. 

It is well known (see e.g., Ref. Is^ that the wave functions 
for free electrons in the reservoir (say, a) with an oscillating 
uniform potential Va (t) are of the Floquet function type: 



(4) 



Here Jn{x) is the Bessel function of the first kind of the n"^ 
order; E = R^fc^/(2me) (me is an electron mass). The cor- 
responding distribution function /o,q = {a[,{E)aa{E)) (here 
(...) means quantum-statistical averaging) is independent of 
the oscillating potential Va and is the Fermi distribution func- 
tion 
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1 + exp 



(5) 



Here Ta is the temperature of the reservoir a; ks is the Boltz- 
man constant. 

In general, to find the Floquet scattering matrix S'p, we 
have to investigate the transmission and reflection amplitudes 
of electrons with a wave function tp{E,t,r) given by Eq.Q 
incident on the oscillating scatterer. However if the frequency 
UJ is small compared with the energy of electrons participating 
in the transport (i.e., with the Fermi energy fi) 



(6) 



we can reduce the problem to scattering of ordinary plain 
waves. To this end we use the following trick^^. We imagine 
that in the leads connecting scatterer to the reservoirs the 
oscillating potentials tend to zero: Va ~ 0. Then in the leads 
the electron wave functions are simply plain waves 



ipo,a{E,r) = 



(7) 



In this region we introduce annihilation operators a, b for in- 
coming and outgoing particles, respectively. In close analogy 
with Eq.Q they are related but through the Floquet scatter- 
ing matrix SF,ai3{E„, E) describing scattering of incident and 
outgoing plane waves: 



6U-B)=^ ^ 5^,„M^.^n)«^(-^5n). (2) 

,a E„>o 

By definition the reservoirs are not affected by the cou- 
pling to the scatterer and thus they are in an equilibrium 



L(-E) = ^^SF,«,3(-B,£„)a^(£;„). (8) 

P n 

Comparing the wave functions Eq.||lJ and Eq.Q we see 
that the annihilation operators a for particles in the leads 
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can be expressed in terms of the annihilation operators a for 
particles in the reservoirs as follows— 



n— — oo 



'a^{E -nhLu). (9) 



The above representation is valid for small frequencies Eq.Q. 
Thus we can put k{E„) ~ k{E) ignoring the terms of order 
Juj/ii and smaller. In other words we ignore the reflection 
at the interface between the region with oscillating potential 
and the region without one. 

Using Eqs.JHJ and @ we calculate the distribution func- 
tions fi°"'\E) = {bi{E)L{E)) for outgoing and fi'^^E) = 
{al{E)a^{E)) for incoming electrons in the leads as follows 



hui 



.fo,aiE ~ nhu), 



(10a) 



= I] J2 S*p^afj{E,Eq)SFMl3{E,Em) 

(3 n,m,q— — QCi 

f^) Jr^+^ f^) e'<'-™'^^/o,/3(S-nM- 



(10b) 



Now the dc current la of spinless electrons, the quantity of 
interest here, flowing from the scatterer through the lead a 
can be expressed in terms of these distributions^^, 



Ia = lj dE [fi°-'\E) - ft\E)] . (11) 



Substituting Egs. lim into Eg. llH we find 



/3 n— — oc 



Huj 



(12) 



Ea. l|12|l is the basic result which allows us to analyze the dc 
currents flowing in the system under consideration. So far we 
put no restrictions on the reservoirs. Different temperatures 
of reservoirs as well as different (stationary) electrochemical 
potentials can by themselves give rise to dc currents. We 
will not consider the most general situation here. Pumping 
in the presence of stationary chemical potential differences is 
investigated by Entin-Wohlman et ali— i— . Here we focus on 
dynamically oscillating potentials. 

In what follows we assume the reservoirs to have equal tem- 
peratures and equal dc-components of electrochemical poten- 
tials but the oscillating reservoir potentials Va can be differ- 
ent: 



trix S{E) changes signiflcantly^. 



III. ADIABATIC APPROXIMATION 

To zero-th order in frequency the elements of the Floquet 
scattering matrix can be approximated by the Fourier coeffi- 
cients So,n of the stationary scattering matrix So, 



So,n{E) 



2tt 



dte'"'^'SoiE,t). 



(14a) 



Ta = To, Mo, a = Mo, a = 1, 



(13) 



In this case the distribution functions entering Eg. (1121 are 
independent of the lead index: fo,a(i3){X) = fo{X), where /o 
is the Fermi distribution function with temperature To and 
chemical potential fio- 

To calculate the Floquet scattering matrix SpiE, E„) one 
needs to solve the time-dependent scattering problem. Gen- 
erally this can be done only numerically (see e.g., Ref i3Qil . 

Here we are interested in the limit of low frequencies. 
In this limit we can use the adiabatic approximation as 
a starting point and can express the Floquet scattering 
matrix in terms of a stationary scattering matrix with 
time-dependent parameters (the frozen scattering matrix): 
So{E,t) = So{E, {P{t)}). Here {P} is a set of parameters 
Pi{t) = Pi, 0-1- Pi, 1 cos(w/:-|-0i), i = 1, 2, . . . , TVp oscillating with 
frequency ui. The scattering matrix 5'o(-B, {P}) describes re- 
flection and transmission of particles with energy E at given 
frozen parameters Pi. This description is valid if the energy 
scale hu! dictated by the modulation frequency is small com- 
pared with the energy scale 5E over which the scattering ma- 



SoiE,t)^ Y e-'""*So,„(S), 



(14b) 



n— — oc 



as foUows^S 



Sf{E„,E) ~ Sf{E,E- 



So,n{E). 



(15) 

However in general this approximation is not sufficient to cal- 
culate the current to order uu. In particular if the oscillating 
potentials Va ^ are applied to the reservoirs then to calcu- 
late the dc current to first order in frequency uj one needs to 
know the Floquet scattering matrix with the same accuracy. 

Note that fortunately in the case of stationary reservoirs 
{Va = 0) there exists a representation (see Eq.JHJ in Ref 30) 
which allows to calculate the dc current (with accuracy of 
order cu) using only the zero order approximation Eq.HlSt . 
In contrast another representation (see Eq.@ in Ref l3(]ri for 
the same dc current requires the knowledge of the Floquet 
scattering matrix with higher accuracy (i.e., to the first order 
in frequency). 

Note that the nonadiabatic corrections to the scattering 
states and the corresponding corrections to the pumped cur- 
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rent were considered in in the limit of a small mod- 

ulating potential. Our approach is valid for an arbitrary os- 
cillating potential since we take into account the effect of all 
the harmonics of the pump frequency lu. 

To calculate the Floquet scattering matrix with an accu- 
racy of order u we generalize the approach used in Ref 1 and 
start from the unitarity conditions for the Floquet scattering 
matris— 



Since Aaa = 0, we can immediately conclude that the reflec- 
tion {a = /3) coefficients are with accuracy of order to defined 
by the first terms on the RHS of Egs. 11711 . This fact justifies 
our representation for the elements of the Floquet scattering 
matrix in Egs. 11711 . 

We next need to determine the off-diagonal elements of A. 
The detailed calculation is given in Appendix^. The central 
result is the relation (valid to first order in tj), 



Sp^a/siEn, E)SF^aj{En, E,n) = 5mo30^, (16a) 



a n — — oo 



Sp^a/siE, E„)SF,'yl3{Em, En) = SmoSct-i- (16b) 

(3 n — — oo 

Taking into account that Egs. 1151 are a zero-th order ap- 
proximation we will seek the first order approximation in the 
following form 



SF{En,E)=So,n[ ^"^^ ] +niuAn{E) + 0{u;^). (17a) 



SF{E,E^n) = So.n { ^^^ " j+;tC^ln(g)+0(c^'). (17b) 

Here An{E) is a matrix of the Fourier coefficients for some 
matrix A[E,t) = A{E, {P{t)}) which is treated as indepen- 
dent of energy on the scale of the order of nhuj; 0{uj^) denotes 
the rest which is at least of second order in frequency lu and 
which we neglect. Note that the first terms in Egs. 11711 should 
be expanded to the first order in lj 



So.; 



E + E±r. 



SoAE)±n.f-^, 



and other terms (of higher order in uj) should be ignored. 

Based on Eo. l2H we will show that Eg. 11711 is, in fact, an 
expansion in powers of hui/SE. Substituting Eos. 1171 into 
Eos. 1161 and keeping the terms of order lo'^ and uj^ we get the 
required relations which can be used to calculate the current 
Eq.jnj. 

In particular the diagonal part (m = 0, /3 = 7) of Egs. 1161 



gives 



E E So%^,„(i5)A„^,n(i5)+c.c. 

a(j3) n— — 00 



T^We E E n\So,cpAE)\''. 



(18) 



n— — 00 

Here c.c. denotes complex conjugate terms. The sign — (+) 
corresponds to the summation over a(/3). 

In what follows we concentrate on the case without mag- 
netic fields and suppose that the stationary scattering matrix 
is symmetric in lead indexes: 

<S'o,Q,3 = 5o,/3a. (19) 

It follows from Eg. 1181 that in this case the matrix A is anti- 
symmetric: 



hiu {sl{E,t)A{E,t) + A'' (E,t)SQ(E,t)) = -r{Sl;So} 



r{sl;So} = in 



dSl dSo 



dt dE 



dSl dSo 
dE dt 



(21a) 



(21b) 



Here i is the imaginary unit. Since the scattering matrix 
is unitary SqSo ~ I (where 7 is a unit matrix) the matrix 

7'{5'o;5'o} is traceless: J2 'PlSl; So}aa = 0. (see Appendix 

a = l 

^for the detailed proof). 

Note Avron et al4fi consider a closely related matrix tl = 
V{So; Sq}. This matrix is the commutator Cl = j-[T,£] of 

the Wigner time-delay matrisS^iS^: f = -ih^Sl, and the 

matrix of the energy shifli— *— >: £ = ih^^Sl- Note however, 
that on the RHS of Eq. I)21a|l the commutator appears with 
a different sequence of 5*^ and S as compared to f2. For this 
reason (and other reasons to be become clear later on, we 
have introduced a separate notation, the Poisson bracket V. 
As we will show (see Eg. 1341 1 the diagonal elements {e/h)Vaa 
of this commutator are just spectral current densities (current 
per energy). 

From Eq.lHTaJ we find an expression for the product of the 
frozen scattering matrix with elements of A, 



AhujRe[So,al3Aal3] 



1 



V{Sl,So}pp-P{Sl;So}c,c. 



(22) 

Note that the scattering matrix So is a unitary matrix of 
dimension Nr. Evidently Eg. 1221 is consistent with Eg. l2Ul . 

Below we use Egs. 1171 and 1221 to evaluate the current 
Eo. p2^ with an accuracy of order u. 



IV. LINEAR RESPONSE ADIABATIC 
CURRENT 

Now we use the adiabatic approximation introduced in the 
previous section and calculate the zero-frequency, dc-current 
Eg. 1121 to linear order in the oscillating potentials T4< — > of 
the reservoirs at finite temperature To- We assume that the 
following conditions hold: 



hcj < fcsTo, (23a) 



Ac 



eVo, < fcflTo. (23b) 

The first inequality {hu <^ ksTo) is relevant for experi- 
ments on adiabatic {uj 0) quantum transport. The second 
inequality defines nothing but the linear response regime. 
In Eg. 1121 the sum over n contains approximately n,nax ~ 
(20) terms. Therefore, huin < eV and because of Eg. l23b|l 
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we have hojn <C ksT. Hence we can expand the Fermi 
function entering Ea. l|12|l . Taking into account Ea. l|18|l this 
expansion (up to second order in uj) is; fo.i3{E — nTiuj) ~ 

fo{E) - nhuj^^§^ + ^n^h'^uj^^^^^. Substituting this dis- 
tribution into Eg. 1121 and take the sum over n. We use the 
summation formulae for the Bessel functions— 

oo 

Jn + m{X)J„ + q(X) = Smq, 

n— — oo 

oo 

Yl nJn+rn.{X)Jn^q{X) — —UlSmq + —{5m(q + l) + 5m{q~\)), 

oo / < 2 \ 

J2 ri^ Jn+m{X)J„ + q{X) — im^ + Smq 
n — — oo ^ ^ 

-X ([m - 0.5](5„(,+i) + [m + 0.5](5™(q_i)) 

+ ^{5m(q+2) +5m(q-1))- 

(24) 



r 

Here the bar denotes the time average X = ^ J dtX{t) 



over a time period T = 22: • — e^/h is the spinless 
conductance quantum; the function 'P{X; Y} is defined in 
Ea. ll2Tbl . To arrive at Ea. ll25ct we used the unitarity condi- 
tion ^ I So^af) I = 1 and the fact that the average potential is 

fi_ '_ 

zero: K>(t) = 0. 

We emphasize that in the above expressions we omitted 
all the terms which are of the second (and higher) order in 
frequency lo and/or in potentials Va- Next we characterize 
briefly the three contributions to the current la- 

The current (j^g to solely the oscillating scatterer. 

It determines the quantum pump effect when the reservoirs 
are stationary. It is the formula obtained by Brouwes— . 

The current I^"'^*^ is a consequence of the rectification 
of ac currents flowing in the system under the influence of 
ac potentials K, applied to the reservoirs. In context of 
pumping this effect was considered by Brouwer in Ref. 
Note that this rectified current depends on the conductances 
Ga/3 = —Go\Sai3\'^ and the Corresponding potential differences 



After that, substituting Egs. 11711 . 12211 and applying the 
inverse Fourier transformation Ea. ljl4b^ . we sum over q and 
m. Finally we represent the dc current la flowing in lead 
a under the action of an oscillating scatterer and oscillating 
reservoir potentials in the following way: 



(25a) 



(25b) 



(25c) 



(25d) 



AVaf3{t) = Vi3{t) — Va{t) iu close analogy with the dc cur- 
rent flowing in response to a dc voltage. We stress that here 
AVai3{t) depends not only on the amplitudes of the corre- 
sponding potentials but also on the phase lag Aipap — ipp — ipa 
as well. In particular if the amplitudes of two oscillating po- 
tentials are equal Va = V/s = Vo then the potential difference 
reads 

AVaf,{t) = -2Vo sin(Av5c/3) sin (tut + ^^fi±^^ . (26) 

This equation [together with Ea. l25cl 1 shows clearly that the 
rectification of ac currents can depend on the phase lag be- 
tween the applied ac potentials and, hence, it can mimic a 
quantum pump effec*~. 

The third term J^'"'-' is novel. Interestingly, as we will see, 
this current renormalizes both jjf^™-^^ and la""^^^ ■ The cur- 
rent 7^'"*' is a consequence of the interference between the 
ac currents produced by the external voltages and the ac cur- 
rents produced by the nonstationary scatterer. Remarkably, 
it is essentially determined by commutator expressions. An 



oo 

ia = JdE (-^^) {e^"-'^ + I'r^'^ + 4""'} , 





j(— )(i5)=^^('M5CM5t(i5,t)) 
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li'"'\E) = E {'P{ShSo}p,3 ~ r{Sl;So}aa - Nrr{Sl^p;So,af3}) ■ 
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" oscillating" scatterer is much richer in physics then expressed 
by Eq. I|25c|l . The expression for is widely used but 

this is only a part of a correct answer. The part is 
due to a rectification of external currents caused by the time- 
dependence of the conductances. The oscillating scatterer is 
much richer: It generates its own ac currents which can in- 
terfere with the external ac currents. This interference effect 
leads to 

Before proceeding we check the current conservation. To 
this end we sum la over the lead index a. Note that each 
of the currents j^"^p^ ^ I^^'^*'^ , and li' is separately con- 
served. This fact supports the current decomposition intro- 
duced above. 

For the pump currents using the Birman-Krein 

formula^^ we find 



{pump) 



Tr 



dSp 
dt 



= ^ln(detSo) =0. 



Here we take into account that the average of a time deriva- 
tive is identically zero: ^^||^ = 0; Tr denotes the trace of a 
matrix: TrS = J2S 



The conservation of the rectification currents X] 



(rect) 



Go ^ (Ka(*) ~ |5o,a/3| = follows from the unitarity 

condition |So,a/3|^ = Y \ So,ai3\^ = 1. 



The current J^'"*^ is conserved as well. Since the matrix 
PIS;!,; So} Ea. ll2TBl is traceless we get from Ea. ll25Ht 



E = ^ E ^/'W E (^{^^ ' ^o}/3/3 - ; So}.. 



2Nr 



I 



To shed more insight onto the nature of the new contribu- 
tion we consider a simple but a quite generic example. 



/('■^^"(i?) = GoT{E,t)[V2{t) - Vi{t)], 



(29b) 



A. Two terminal single channel scatterer 

Consider a nonstationary scatterer connected to only two 
reservoirs q = 1,2 via single channel leads. For such a scat- 
terer, assuming there are no magnetic fields, the stationary 
scattering matrix So is a symmetric 2x2 unitary matrix. 



So = e'- 



iVT 



(27) 



Here R and T are the reflection and the transmission prob- 
ability, respectively {R + T — 1). The phase characterizes 
the asymmetry between the reflection to the left and to the 
right. The phase 7 relates to the change of the overall charge 
SQ on the scatterer (for instance a dot) via the Friedel sum 
rule^*: ^7 — ivSQ/e (where e is the electron charge), or in 
different notation SQ — e/(27ri)(5[ln det S]. We assume that 
R,T = 1 — /?, 7 are functions of the electron energy E 
and the external parameters Pi{t) varying with frequency ui. 
Before proceeding we remark that for the case Nr = 2 the 
current J^'"*' Eq. <25dll can be simplified 
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(28) 



-K.(t)P{So%,;So,..} , 



Substituting the scattering matrix Eq.JJZJ into Eas. H25^ and 
1281 1 we find the currents 7i = —I2 fiowing between the scat- 
terer and the reservoirs: 



{pump) 



(E)^—R{E,t) 



de(E,t) 

dt ' 



(29a) 



(29c) 



4n \ at BE 



lit) [^2(t) 



Vi{t)]. 



These expressions demonstrate that the current has 
common features with both the rectification current /(''^'^*) 
and the pumped current j(p"™-p\ Like the former, the cur- 
rent depends on the potential difference AV12. Like the 
latter, the current can exist even at equal reservoir po- 
tentials Vi(t) = V2{t). In this case, the conditions necessary 
for the existence of j'^"*) and /(f"™") are the same^: First, 
the scatterer has to be asymmetric, i.e., 7^ 0, and, second, 
the time reversal symmetry (TRS) has to be broken. We note 
that the current depends on both the oscillating reser- 
voir potentials Va{t) = Va cos{ujt + ipa) and the oscillating 
scattering parameters Pi{t) — Pifi -f Pi^\ cos{ut + There- 
fore analyzing the presence/absence of the TRS we have to 
consider all the phases, namely as well as 4>i- 

We have here treated only non-interacting electrons. As a 
consequence sums of potentials appear in Ea. ll29cl . This is 
in contrast with an electrically self-consistent theory which 
permits only the appearance of voltage differences. If inter- 
actions are switched on- then the (self-consistent) potential 
[/ 7^ inside the scatterer becomes dependent on external 
potentials and the differences Va — U should appear in- 
stead of v.. U is in general a function of all the oscillating 
parameters Pi(t), all the external potentials Va and also of 
the potentials at the gates which influence the electrostatic 
potential inside the scatterer. Our expressions do, however, 
conserve current. 

We see that the flrst term on the RHS of Ea. l29cll renor- 
malizes the pumped current /(p"™?) and the second one renor- 
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malizes the rectification current 7(^'=='). The latter is due to 
nonadiabatic (first order in uj) corrections to the conductances 
arising from the corresponding corrections Eq. 11711 to the scat- 
tering matrix._Note that the analogous corrections are dis- 
cussed in Refs ;2S:46 in context of pumping in the presence of 
a dc bias. 

Since the pump effect is the main topic of this work we 
consider now the case with Vi{t) = V2{t) in more detail. 
This case corresponds to an experimental setup in which the 
scatterer and a large portion of the reservoirs to which it 
is connected are subject to long wavelength radiation. The 
effect of such radiation can be modeled via an oscillating uni- 
form potential V{t) which is the same at different reservoirs: 
Vi{t) — V2{t) = V{t). In this case the rectification current 
Ea. l29bl is absent = 0, and the whole dc current la 

can be reduced to the simple form 



afo(E) 

SB 



R{e,t) 



de{£,t) 
dt ' 



(30) 



£ = E + eV(t). 



To obtain this result we have used the following identity: 
-A|S = rSA_ ^-^^Yi A = ey||. We have also introduced 
the full time derivative: 4, = 4: + e^-^- 

dt at dt at 

This result can be understood in the following way: For 
stationary reservoirs {V{t) = 0) the pumped current is de- 
scribed by equations 125all and Il29al with the quantities R 
and 6 taken at the energy E of incident electrons. However, 
if the chemical potential fi{t) = ^o + eV(t), V(t) 7^ oscillates 
slowly (u ~* 0) then we can consider incident electrons having 
energy £ = E -\- eV(t) following adiabatically the reservoir's 
potential V{t). Substituting in Ea. H29a|l £ instead of E and 
replacing a partial time derivative by a full time derivative we 
get Eg.CT. 

ft should be noted that the above substitution £ — E + 
eV{t) implies that the potential inside the scatterer {U = 0) 
is independent of the external potentials Va . This is correct 
for noninteracting electrons but it should be modified if the 
interactions are present^. 

From equation 13UII we can conclude that the effect of an 
oscillating external potential V{t) is like the effect generated 
by an oscillating parameter of the scatterer (i.e., an oscillat- 
ing internal potential) . Therefore to analyze the ability of an 
open system (the scatterer plus reservoirs) to generate adia- 
batic dc currents we have to consider the full set of oscillating 
parameters {Vc{t),Pi{t)},{a = 1,2, . . . , Nr;i = l,2,...,Np). 



V. DC CURRENT AT FINITE AC VOLTAGES 

Now we go beyond linear response theory. We suppose that 
the potentials Va can be large compared to the temperature. 
Thus we calculate the current Eq. 1121 with accuracy up to the 
first order in uj and with an arbitrary ratio of the potentials 
Va to the temperature: 



eVa < /iO.a- 



(31a) 



(31b) 



Since the potentials Va are not necessarily small compared 
with the temperature T we can not expand the Fermi function 
fo./iiE — nhui) entering Eg. 11121 . Nevertheless Ea. ll31al allows 
us to sum over n and to simplify Eg. 11211 . 

To this end we go from the energy representation over to the 
time representation. We express the Fermi function fo,i3{E) 
Eq.Q and the Bessel functions Jn{x) as follows: 



00 

foAE-nhco)^ J dr/o,^(r)e'(^-"''"'S, 



i 

^ J dtW^{t)e- 



Jn-\ 



eVfi 



T 



dtiVP^^(ti)e'<"+'"''"*i, 



foAr) = 



2;isinh (TrfcsTo^) 



-'1^0. 



Substituting these equations into Eg. 1)12^ and summing 
over n we obtain a delta-function 5{ti —t — t) which allows us 
to perform one additional integration. At r > (r < 0) we 
integrate over ti (t). This leads to the substitution ti = t + r 
{t = ti + \t\). Further we expand sin(wT + uit + (pp) to first 
order in ut. We can do this because for the relevant r < 



Eq. 131 at gives lot <^ 1. Next we integrate over r and finally 
get the dc current as follows: 



Ia^jJ'dE^Jdtlj2 E fo{E + {q + rn)^-,Mi))shc,0{E,E,)Sp,a,3{E,E^)e'^'--''^-' -fo{E;^^^^^ (32) 



/3 m,g — — 00 



Here we have introduced the Fermi function with time- 
dependent chemical potential fia{t) = /io.a -I- eVait) Eq.l|^: 

/o(S;M.(i)) = [l + exp(^^)] \ 

Note that Eg. 13211 is valid both for the adiabatic as well 
as for the nonadiabatic case. The only restriction is that the 
frequency has to be small compared with the temperature 



Eq.lEIaJ. 

Next we use the adiabatic approximation of Sec lllll and cal- 
culate the current la to first order in frequency uo under the 
conditions of Eg. 11311 . To this end we substitute Egs. 11711 and 
I22II into Eg. 1321 1 and expand the Fermi function in powers of 
UJ. Next we use the inverse Fourier transformation Eg. ll4bll 
and after a little manipulation (we integrate by parts on en- 
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ergy and dropped the contribution arising from _B = 0; in 
addition we exploit the unitarity of the frozen scattering ma- 



trix ^ I So,a/3(-B, t)| = 1) and find the current: 



J 



2Nr 



(33) 



The above equation generalizes Eqs.JJSJ to the case of finite 
voltages. Current conservation Ic, — can easily be proven 

a 

in analogy with Eas. l)25^ . 

Next we concentrate on the pump effect and consider 
the case with reservoirs having equal oscillating potentials: 
lic,{t) = nit) — fj.0 + eV co8{ujt + tp), a = 1, . . . ,Nr. Since 
the Fermi functions entering Ea.l|88|l become independent of 
the lead index we can sum up over f5 and obtain 



tx> r 
/c = J dE^ Jdtfo{E;ii{t)) 
b 



dIa(E,t) 
dE ■ 



(34) 



dE ^ h' \"^0,»3oJ<»« — '27r I dt dE dE dt 



Here we used an obvious equality: 'P{SQ;So}aa ~ 

The quantity dia {E, t) / dE is the spectral current density 
at energy E and time t, (i.e., the current within the energy 
interval dE) produced by the adiabatically evolving scatterer 
towards the reservoir a. This definition seems reasonable be- 
cause of a conservation law ^ dIa/dE{E, t) = which is valid 

at any energy E and at any time moment t. Note that in the 
case of stationary reservoirs the same interpretation was given 
in Ref. li^. 

These currents (or more precisely, the ability to produce 
them) are an intrinsic property of a time-dependent scatterer. 
This property differentiates between a nonstationary scatterer 
and a "frozen" one. Note that the Fermi distribution func- 
tion in Eg. (1341 describe the filling of (potentially) existing 
"current" states of a nonstationary scatterer. 

At V = the Eg. 13411 reproduces Brouwer's result Eg. ll25bll 
and agrees with that obtained in Refill. At small voltages 

— > for the scattering matrix Eg. l27ll we get Eg. l3UII . 

Equation 1341 determines the dc-current to the first order 
in u> pumped by the slowly oscillating scatterer between the 
reservoirs having equal (possibly zero) oscillating potentials 
Va{t) — V{t). Formally in the adiabatic case under con- 
sideration the effect of oscillating chemical potentials is only 
the change of an energy of electrons falling upon the scat- 
terer. However, in fact, the phase ip of an oscillating potential 
V{t) — V cos{ijjt + tp) is of a great importance because of the 
following. An adiabatically pumped current Jc, 7^ is gener- 
ated already if the time reversal symmetry is broken in the 
whole system including the scatterer and the reservoirs. At 
V analyzing this question we have to take into account a 
possible phase shift between the potentials of reservoirs and 
the oscillating parameters Pi{t) — Pi,o + Pi^i cos{Lut + (j>i) of 
a scatterer. In particular even a scatterer with a single oscil- 
lating parameter can produce an adiabatic dc current if only 



VI. INSTANTANEOUS CURRENT 

In this Section we derive an expression for the instan- 
taneous current of an adiabatic quantum pump simultane- 
ously subject to oscillating external potentials. We first clar- 
ify the physical meaning of the (diagonal elements of the) 
quantity 7'{S'q;S'o} defining (antisymmetric in lead indices) 
nonadiabatic corrections Egs. 12111 to the scattering matrix. 
From the geometrical point of view" this quantity is a cur- 
vature in the time-energy plane. The physical interpreta- 
tion is based on Eg. 1341 . We can consider the quantity 
dIa/dE{E,t) — j^V{So; SQ}cia as an instantaneous spectral 
current which is pushed by the oscillating scatterer into the 
lead Q. A more detailed partitioning of the current follows 
from Eg. 1331 : We can say that the scatterer drives the fol- 
lowing spectral currents from the lead /? into the lead a: 

e r{So; Sl}cc. " r{So; S^jp^ + iV.P{go.c^; gp^^} 



dl, 



af3 



dE 



2Nr 



(35) 

The above spectral currents are subject to the following con- 
servation law: X] dI{E,t)ap/dE — 0. This property sup- 

a = l 

ports the point of view that these currents arise "inside" the 
scatterer (they are generated by the nonstationary scatterer) 
without any external current source. Thus we can consider 
the pump as a source of currents rather then a source of 
voltages^'^^ 

For a scatterer with scattering matrix Eg. 1271 . we obtain 
the spectral currents 



dIii{E,t) 
dE 



dl22{E,t) 

dE 



e_ / a(7-6>) dR _ djj - 9) dR 
i-Ky dt dE dE dt 



e_ f d{-f + e) dR 
'47r I dt dE 



9(7 + 6*) dR 
dE dt 



(36a) 



(36b) 



Note that above currents depend on the phase 7 related to 
the charge of a scatterer. 

Strictly speaking if we are dealing with time dependent 
currents (instead of only the time averaged currents) then 
we need to show that these currents satisfy the continuity 
equation for the charge currents: 



(37) 



Here Ia{t) is the full time-dependent current fiowing through 
the scatterer to the lead a; Q{t) is a charge of a scatterer. 
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To calculate Ic, (t) we first calculate the Fourier transformed 

T 

current I^.i = ^ / dte''-"* Ia{t) which reads^: 



dE 



\^{bi{E)bc{E + Ihu;)} - {ai{E)ac{E + Ihij)}^ 



The operators ba and aa are defined in Eqs.(|HJ and @, re- 
spectively. The calculations analogous to those leading to 
Eg. 1331 give us Ia,i. Performing the inverse Fourier trans- 
formation Ea. lll4bll we finally get the time-dependent current 
Ia{t) flowing in the system as follows: 



Ut)^ I dEj2\l[fo{E;Mt))-foiE;Mt))]\SoME,tf-e-^^ 



h[E-iip{t)) 



dN^p{E,t) 



dE 



+ h{E;t.,{t)f-^§^ 



Here we have introduced the partial density of states for a 
"frozen" scatterer, 

dNap _ i ( dSl^afj dSo,al3 \ 

~dE~ ~ 4^ Iv^^ " °'''^^^) ■ 

These density of states define the charge Q{t) of a "frozen" 
scatterer as follows: 

oo 

QW-^EE / dEf,{E;^.,{t)) ''''Zf''^ (39) 

-/So 

The quantity Ia{t) Eq.l^HJ and Q{t) Eq.lEHll do satisfy the 
continuity equation 1137^ . 

The three terms in the curly brackets on the RHS of Ea. 13811 
can be interpreted as follows. The first term defines the cur- 
rents flowing under the action of external voltages through a 
"frozen" scatterer. The second one defines currents attributed 
to an oscillating charge of a "frozen" scatterer. The third term 
can not be entirely viewed just as a nonadiabatic correction 
of either the "frozen" conductance nor of the "frozen" density 
of states. It is more naturally to consider it as the ac currents 
generated by the oscillating scatterer. The ability to gener- 
ate these ac currents differentiates a nonstationary dynamical 
scatterer from a merely "frozen" scatterer. 



VII. DISCUSSION 



the oscillating potentials Va it) of reservoirs can be accounted 
for by allowing the energy E of incident particles to follow 
adiabatically the reservoir potential: _E — > f = _B + eVait). 

We emphasize the importance of the phases of all the cycli- 
cally evolving quantities (the potentials of reservoirs and the 
parameters of a scatterer) for generating a dc current. In par- 
ticular even when all the reservoirs have the same oscillating 
potential Va{t) = V{t) and the rectification effect is ineffec- 
tive the dc currents &t V = Q and at 1/ 7^ can nevertheless 
differ significantly. 

The analysis allows us to perform a current partition that 
clarifies the physical meaning of the (diagonal elements of the) 
quantity ■P{S'q;S'o} (the Floquet corrections to the "frozen 
scattering" matrix) and show that they correspond to spec- 
tral current densities generated by a dynamic scatterer. The 
instantaneous current contains a contribution from such self- 
generated ac currents in addition to the currents from the 
"frozen" charge and the ac currents generated by the exter- 
nal potentials. 

We emphasize that the results presented in this work, the 
effect of external ac-potentials on a quantum pump, are of 
importance whenever the pump is not part of an ideal zero- 
impedance external circuit. In particular, if the pump is in 
series with a resistance used to measure the voltage generated 
by the pump, or if the circuit is a multiterminal circuit with 
probes used to measure voltages, the results presented here 
will be needed. 



We have investigated the nonstationary adiabatic charge 
transport through a time-dependent mesoscopic scatterer cou- 
pled to reservoirs subject to oscillating voltages. The external 
voltages applied to the reservoirs induce ac currents flowing 
through the scatterer. In addition the oscillating scatterer it- 
self is a source of ac currents flowing between the reservoirs. 
In general these two types of currents interfere with them- 
selves. This gives rise to renormalization of the rectification 
(i.e., proportional to the potential difference) contribution to 
the dc current and gives rise to a renormalization of the quan- 
tum pump current. 

To analyze this interference effect we calculated the Floquet 
scattering matrix beyond the adiabatic approximation. We 
investigated the first order in uj corrections to the (adiabatic) 
scattering matrix and found that the dc currents of both the 
zeroth and the first order in lu can be expressed in terms of a 
stationary scattering matrix with time-dependent parameters. 
Within this approximation, within a non-interacting theory. 
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APPENDIX 
1. The matrix A 

The matrix A defines the first order in frequency corrections 
to the adiabatic Floquet scattering matrix Eos. 1171 . 12U1 . It 
is a matrix antisymmetric in lead indexes. 

To obtain Eg. 1211 we substitute the adiabatic expan- 
sion Eg. ll7al in the current conservation condition Eg. ll6al . 
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Keeping terms of order oj" and lu^ we get the following: 

S*F,afl{En, E)Sp,a-i{En, Em) 



And second, we can write q%q^ 



ingly, 



['S'd'S'o 



and, correspond- 



ed n — — oo 



' 9£ 9t dt dE 



X S, 



n + m dSo. 



(E) 



2 dE 

E E ^0,aP,n{E)So,ai,n-m{E) 
a n 

n + m dSo,a',,n-miE) 



+ hujA 



a-y,n~7n 



(E) 



Substituting Eas. lA.21 and IIA.311 in Eg. llA.lll we arrive at the 
Eq.lHTJ. 



U,Q:'y,ri — m V 

dE 



+ E E ( n.^^^^^^ + nu:A*^,,„iE) ] So 



+ TLUjAaj,n-m{E) 

(E) 



Applying the inverse Fourier transformation Eq. l|14b|l and in- 
troducing corresponding matrixes we rewrite above equation 
as follows: 

(^\So{E,t)\^ + hujSl{E,t)A{E,t) + hu;A\E,t)So{E,t) 



2. The commutator matrix V 

The matrix ViS^; So} defined in Ea. ll21bll is self adjoint 
r{Sl;So}=V^{Sl,So}, (A.4) 

and traceless 



Tr[P{Sl;So}] =0. 



(A.5) 



_L « c-t a'' So _ S. r At SSo 1 

2 "^0 dEdt '"et ['^0 dE J 2 asat 



iS-^So 



To demonstrate the latter property we use the equality 
d[S']S'^ = —Sd[S^] following from the unitarity of the scatter- 
ing matrix SS^ — I and the invariance of trace to the cyclic 
rearrangements of the matrices. As a result from Ea. ll21b| we 
get 



— m 



Tr[V] = ihTr 



(A.l) 

To simplify further this equation we use the unitarity condi- 
tion for the frozen scattering matrix: So{E,t)SQ{E,t) = /. 
First, from this condition it follows that 



at dE as ' 



at 



So{E,t)\' 



/3"f, — m 



iA^ 



ihTr 



= ihTr 



sslast^ _ et aso ssl a 
at as ^ as at 



at as at aE 



= 0. 
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